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Double parton scatterings in high energy hadronic colhsions 

G. Calucci and D. Treleani''* 
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Trieste, I 34014 Italy 

CDF has recently measured a large number of double parton scatterings. The observed value of (Jeff, the 
non perturbative parameter which characterizes the process, is considerably smaller as compared with the naive 
expectation. The small value of (Jeff is likely to be an indication of the importance of the two-body parton 
correlations in the many-body parton distributions of the proton. 



1. INTRODUCTION 

The inclusive cross section for a double par- 
ton scattering, namely of an event where, in the 
same inelastic interaction, two different pairs of 
partons scatter independently with large momen- 
tum transfer, is written aslOll: 



Pt 



D2{xA,x'^;h)a{xA,XB) 

a{x'^, x'g)D2ixB , x'g;h) 
dhdx Adx sdx'j^dx'g 



(1) 



ct{xa,xb) is the parton-parton cross section in- 
tegrated with the cut off p^, which is the lower 
threshold to observe final state partons as mini- 
jets, X is the momentum fraction, A and B are 
labels to identify the two interacting hadrons. ub 
is a function of the product d{xA,XB)<j{x'A,x'g). 
Actually the two different partonic interactions 
are localized in two regions in transverse space 
with a size of order {l/pff' and at a relative dis- 
tance of the order of the hadronic radius r, in such 
a way that the two partonic interactions add in- 
coherently in the double scattering cross section. 
The non perturbative input in Eq.(l) is the two- 
body parton distribution D2{x,x' ]\)), which de- 
pends on the fractional momenta of the two par- 
tons taking part to the interaction and on their 
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relative transverse distance b. The transverse dis- 
tance b has to be the same for the two partons 
of hadron A and the two partons of hadron _B, in 
order to have the alinement which is needed for 
a double collision to occur. D2 is a dimensional 
quantity and therefore the process introduces a 
non perturbative scale factor which is related to 
the hadronic transverse size. 
The simplest possibility to consider is the one 
where the dependence of D2 on the different vari- 
ables is factorized: 



D2{x,x'-h)^feff{x)feff{x')F{h) 



(2) 



Jeff is the effective parton distribution, namely 
the gluon plus 4/9 of the quark and anti-quark 
distributions and F(h) is normalized to one. Mul- 
tiparton distribution are then uncorrelated and 
D2 does not contain further informations with re- 
spect to the one-body parton distribution (actu- 
ally feff) apart form the dependence on b, whose 
origin is the dimensionality of D2 and which gives 
rise to the scale factor Ueff- In fact in this case 
one may write 



with 
1 

and 



^eff 



J F^{h)dH 



(3) 



(4) 



feff{xA)feff{xB)^{xA,XB)dxAdxB, (5) 
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the single scattering expression of the perturba- 
tive QCD parton model. 

Eq.(2) is the basic hypothesis underlying the 
signature of a double parton collision which one 
has been looking for experimentally ||J^. The ex- 
pected characteristic feature of a double collision 
is in fact that it should produce a final state anal- 
ogous to the one obtained by super-posing two 
single scattering processes. By looking at the de- 
pendence of CTe/ / on a: CDF has been able to ver- 
ify the correctness of the factorization hypothesis 
in Eq.(2). The range of values of x available is 
limited to x < .2, for the interaction producing 
a pair of minijets, and to x < .4 for the inter- 
action giving rise to a minijet and a photon. In 
the limited range of values of x available, the fac- 
torization hypothesis has shown to be consistent 
with the experimental evidence. 

Since the uncorrelation hypothesis does not 
contradict the experiment, one can work out the 
case where all multiparton distributions are un- 
correlated and one may look for the sum of all 
multiparton interactions to the hadronic inelas- 
tic cross section. The subset where all multiple 
parton collisions are disconnected can be easily 
summed up in the uncorrelated caseQ. The re- 
sult is the semi-hard hadronic cross section an, 
which represents the contribution to the hadronic 
inelastic cross section from events with at least 
one semi-hard partonic interaction. The actual 
expression is 



(6) 



The integration on the impact parameter of the 
hadronic collision, gives the dimensionality to 
the cross section. The argument of the integral 
has the meaning of a Poissonian distribution of 
multiple semi-hard partonic interactions with av- 
erage number depending on the impact parame- 
ter. 

The actual value of Ce// can be obtained by 
taking twice the opposite of the second term of 
the expansion of an in powers of multiple colli- 



sions, so the actual value of an is related to the 
value of (Te// through Eq.(6). The single and the 
double parton scattering cross sections are how- 
ever related to the average number of parton scat- 
terings and to the second moment. Indeed if one 
writes the average number of parton scatterings 
one obtains: 



(n)(JH 



^-^ TV. 

n=l 

J d^pasFip) = as (7) 



and for the second moment: 



{n{n - 1)) aH 



(n(n-l)a5f(/?))\ _.,^(^) 



d'p al [Fipf^ao 



(8) 



The relation between as and (n) and the relation 
between ao and {n{n — 1)) just obtained do not 
hold only in the simplest case of the Poissonian 
distribution of multiple parton collisions. They 
are indeed much more general validity. One can 
in fact obtain the same relations in the most gen- 
eral case of multiparton distributions and keep- 
ing moreover into account all semi-hard parton 
rescatteringsH. One may therefore write 



{n)aH = as and {n{n — l))aH = an 



(9) 



The effective cross section is defined by the rela- 
tion 



(10) 



ao = 

one may therefore write 
\2 o-ff 



{n{n - 1)) = {ny 



^eff 



(11) 



which implies that in case of an overall Poisso- 
nian distribution of multiple parton collisions one 
would have fie// = o'h- When the number of par- 
ton collisions is very large, in the simplest uncor- 
related case, the distribution is Poissonian at a 
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fixed value of tlie impact parameter. Tlie expec- 
tation is tlierefore tliat the overall distribution 

in the number of parton collisions has a larger 
dispersion as compared with the Poissonian case. 
In that regime Ue// is therefore smaller with re- 
spect to an- The comparison between the actual 
value of Geff and of gh depends on the func- 
tional form of F{P). In the simplest case where 
F{l3) = exp(-/3^/i?,^)/7ri?^ one obtains a closed 
analytic expression for au'- 

CTff = 27ri?,2[7 + lnK-|-£;i(K)] (12) 

where 7 = 0.5772 ... is Euler's constant, k = 
as/i^TrR^) and Ei{x) is the exponential integral. 
In this example the relation with the hadronic 
radius r is i? = ry/2. For small k one ob- 
tains gh — * = as, for large k, namely 
tJs — > 00, one obtains (Th 2iTR^{"f + Inti). 
Here aeff = 2-1: E?. The value of an is there- 
fore proportional to the measured value of cTe//, 
the proportionality factor is slightly dependent 
on energy and on the cutoff. Sensible values of 
the hadron-hadron cm. energy and of the cut- 
off give values for an which are some 30 — 40% 
larger with respect to the value of a^f / . Different 
analytic forms for F{(i) give qualitatively similar 
results. 

The effective cross section quoted by CDF is in- 
deed different with respect to the effective cross 
section discussed here and in most of the papers 
on double parton scatterings. (Jeff has a simple 
link with the overlap of matter distribution in the 
hadronic collision when (Jd is obtained from the 
second moment of the distribution in the num- 
ber of partonic collisions, as discussed above. In 
the sample of events with double parton colli- 
sions CDF on the contrary has removed all events 
where triple parton collisions are present. The 
correction is not a minor one since the fraction of 
events with triple collisions is 17%. In the sim- 
plest uncorrelated case just discussed the double 
parton scattering cross section measured by CDF 
would correspond to the expression 

[Mcn. = jd'^^^^^^-'''''' (13) 

The relation above shows which is the complica- 
tion introduced by the requirement of an exclu- 



sive cross section. In order to make the exclu- 
sive selection of the events with double parton 
collisions only, one has to introduce the exponen- 
tial factor which represents the probability of not 
having any further parton interaction. This fac- 
tor, in principle, has a rather complicated depen- 
dence on the overlap of the matter distribution 
of the two hadrons since, by unitarity, it is re- 
lated to the whole series of multiple parton colli- 
sions. The effective cross section quoted by CDF, 
{<Jeff)cDF = 14.5 ± 1.7t2'3™^> refers to the ex- 
clusive measurement and therefore it has to be 
regarded as an upper bound on the value of the 
effective cross section related to an inclusive mea- 
surement, as it has been presently discussed. 

The experimental indication is therefore that 
the effective cross section is rather small as com- 
pared with the naive expectation. The simplest 
assumptions underlying the derivation of Eq.(6) 
have therefore to be revised. 

The main hypothesis which has been done to 
obtain the expression for an in Eq.(6) is the Pois- 
sonian multiparton distribution. On the other 
hand one has to expect correlations between par- 
tons as a consequence of the binding force. While 
most probably correlations will affect the x de- 
pendence of the multiparton distribution only for 
finite values of a;, and therefore at large rapidities, 
correlations in the transverse parton coordinates 
are present in every kinematical regime. Indeed 
the main reason of interest in multiple parton 
collisions, besides the identification of the pro- 
cess itself, is precisely the measure of the many- 
body parton correlations, which is an information 
on the hadron structure independent on the one- 
body parton distributions usually considered in 
hard processes. 

In the next paragraph we discuss the most gen- 
eral expression for the semihard cross section cth, 
which one obtains by 

1) assuming that only two-body parton correla- 
tions are present in the many-body parton distri- 
butions and by 

2) summing all disconnected multiple parton in- 
teractions. 
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2. SEMI-HARD CROSS SECTION AND 
CORRELATIONS 

At a given resolution, provided by the cut off 
p-nim ^^^^ defines the lower threshold for the pro- 
duction of minijets, one can find the hadron in 
various partonic configurations. The probabil- 
ity of an exclusive n-parton distribution, namely 
the probability to find the hadron in a configura- 
tion with n-partons, is denoted by Wn{ui . . . Un)- 
Ui = {hi.Xi) represents the transverse partonic 
coordinate and longitudinal fractional momen- 
tum Xi while color and flavor variables are not 
considered explicitly. The distributions are sym- 
metric in the variables Ui. One defines the gener- 
ating functional of the multiparton distributions 
as: 



J{ui) . . .J{Un) 



Wn {Ui . . .Un)dui 



. du„ 



(14) 



where the dependence on the infrared cutoff p™*" 
is implicitly understood, and one may introduce 
also the logarithm of the generating functional: 
T[J] = ln(Z[J]). The conservation of the proba- 
bility yields the overall normalization condition 



Z[l] = 1. 



(15) 



One may use the generating functional to derive 
the many body densities, i.e. the inclusive distri- 
butions Dn{ui . . . Un)'- 



Di{u) 

D2(ui, U2) 



5Z 

5J{u) 



j=i 



5J{ui)5J{u2) 



j=i 



(16) 



The many body parton correlations are defined 
by expanding J-[J] in the vicinity of J = 1: 



T[J] 



D{u)[J{u) - l]du 



+ / C„(wi...u„)[j(ui)-l] ... 

n=2 ■ 

... [j{Un) - l]dui . . .du„ (17) 

Here D — Di and the correlations C„ describe 
how much the distribution deviates from a Pois- 
son distribution, which corresponds in fact to 
C„ = 0, 71 > 2. 

In the case of hadron-nucleus and nucleus- 
nucleus collisions a systematic use of the AGK 
cutting rules [|j allows one to express the total 
inelastic cross section as a probabilistic super- 
position of nucleon-nucleon interaction probabil- 
ities!^. The same feature holds for the self- 
shadowing cross sections]^. When considering 
hadron-hadron collisions as interactions between 
objects composed with partons, one can make the 
assumption that similar relations hold with nu- 
cleons in the place of nuclei and partons replac- 
ing nucleons. Of course, contrary to the nucleon 
number in the nucleus the parton number is not 
fixed. In this respect semihard parton-parton in- 
teractions have to be regarded as a particular case 
of self-shadowing interactions |9|. The semi-hard 
nucleon-nucleon cross section is then expressed as 
the sum of all the probabilities of multiple parton 
collisions: 



d^PaniP) 



(18) 



with 



/ n! oJiui 

n ^ 
y-^ _! 6_ 



) 6J{Un) 

6 



Za[J] 
Zb[J'] 



m\ 5J'{u[ -P)'" SJ'{u'^ - (3) 

71 'irt 

{ 1 - n n [1 - '^'^'1 («' } n 



1=1 3=1 



.J=J'=0 

(19) 



where /3 is the impact parameter between the two 
interacting hadrons A and B and (Ti.j is the ele- 
mentary probability for parton i (of A) to have a 
hard interaction with parton j (of B). The semi- 
hard cross section is constructed summing over all 
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possible partonic configurations of the two inter- 
acting hadrons (the sums over n and m) and, for 
each configuration with n partons from A and m 
partons from B, summing over all possible mul- 
tiple partonic interactions. This last sum is con- 
structed asking for the probability of no inter- 
action between the two configurations ( actually 
nr^illjlill ~ ^ij] )■ One minus the probabil- 
ity of no interaction is equal to the sum over all 
semi- hard interaction probabilities. 

The presence of multiple parton interactions is 
induced by the large flux of partons which is effec- 
tive at large energies. The most important con- 
tribution to the semi-hard cross section, as a con- 
sequence, is the contribution of the disconnected 
partonic collisions, namely the interactions where 
each parton undergoes at most one semi-hard col- 
lision. These are, in fact, those multiple partonic 
interactions that, at a given number of partonic 
collisions, maximize the parton flux. Indeed the 
search and the observation of the first evidence of 
multiple semi-hard parton interactions has been 
focused to the case of double disconnected par- 
ton interactions H,^. We simplify therefore the 
problem by expanding the interaction probabil- 
ity ( the factor in curly brackets ) as sums and 
by removing all the addenda containing repeated 
indices: 



71, m 



(20) 
(21) 



as a result the semi-hard cross section is con- 
structed with multiple disconnected parton col- 
lisions only, where disconnected refers to the per- 
turbative component of the interaction. Because 
of the symmetry of the derivative operators in 
Eq.(19) one can replace the expression in Eq.(21) 
with: 

nmCTii — -^n'{n — l)m(m — 1)o'ii(T22 + . . . (22) 

in such a way that the sums over m and n can be 
performed explicitly. As a consequence the cross 



section at fixed impact parameter, aniP), can be 
expressed by the operatorial form: 



aH{(3) = 1- exp{-S ■ a ■ S') 
■ Za[J + 1]Zb[J' + 1] 



J=.J'=0 



(23) 



We have avoided writing explicitly the variables 
u and u' and the functional derivative 5/5J{ui) 
has been simply indicated as 5i. 

The form of (Jh{P) given by Eq.(23) is still 
too complicated to be worked out in its general 
form, since all possible multi-parton correlations 
are present in Z. Therefore we further simplify 
the problem by taking into account two-body par- 
ton correlations only. Our explicit expression for 
T is therefore: 



•^[•^+1]= / D{u)J{u)du 

+ i / C{u,v)J{u)J{v)dudv (24) 

where D{u) is the average number of partons and 
C(m, v) is the two-body parton correlation. 

Either by using techniques of functional inte- 
gration or by means of a suitable diagrammatic 
expansion jl^ one is able to obtain in this case a 
closed expression for aniP)'- 



an (/?)=: 1 - exp 



(25) 



where a„ and &„ are functions of the impact pa- 
rameter P and are given by 



DA{ui)a{ui,u[)x 

X Cb{u'i - /3, M2 - /3)(t(u2,U2)C'a(u2,U3) • ■ ■ 

■■■ Deiu'^- P)Y[duM (26) 



bn = 



Ca(w„,ui)<5-(ui,u'i)x 
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X Cb{u'i- (3,U2- f3)a{u'2,U2)... 

■■■ C's«-i /?)a-«,«n) 

n duidu'i . (27) 

The actual expression for a„ holds for n odd. 

When n is odd one may also have the symmetric 
case, where the expression begins with Db and 
ends with Da- When n is even the initial and 
final distribution are either both D.y or both Db- 
In the definition of &„ n is always even, so that 
one of the ends is A and the other is B- One 
may notice that, at a given order in the number 
of partonic interactions, one can obtain a term of 
kind a from a term of kind b by replacing one C 
with a pair of _D's. The operation can be done in 
n ways. The combinatorial meaning of the 1/n 
factor multiplying each term of kind b in Eq.(25) 
is then miderstood. The factor 1/2 in Eq.(25) is 
the consequence of the symmetry between A and 
B. 

The cross section is given by an integral on the 
impact parameter of the interaction probability, 
(Th(/3), that is expressed as one minus the prob- 
ability of no interaction. The probability of no 
interaction is given by the negative exponential 
of the sum over all possible different connected 
structures, namely all structures of kind a„ and of 
kind bn- With our approximations, Eq.(21) and 
Eq.(24), these are in fact all possible connected 
structures which can be built with the average 
numbers Da,b, the two-body correlations Ca,b 
and the interaction a . Expanding the exponen- 
tial, the cross section can then be expressed as the 
sum over all possible structures, both connected 
and disconnected. 

One will notice that, when no correlations are 
present, all terms of kind b disappear and only 
the first of the terms of kind a, namely Da^Db 
is left. In that limit the cross section is given 
simply by: 

aH=yd'/3{l-e-<"W>} (28) 
where 

(n(/3)) = [ Da{u- P)DB{u')a{u,u')dudu' (29) 



which corresponds to the Poissonian distribution 
discussed in the introduction. 

3. TWO DIFFERENT QUALITATIVE 
FEATURES OF THE CORRELATION 
TERM 

The small value of (Teff, the dimensional pa- 
rameter characterizing double parton scatterings, 
which has been measured recently by CDF, is an 
indication that two-body parton correlations, in 
the many-body parton distribution of the proton, 
are likely to be sizable. In the case of an uncor- 
related many-body parton distribution, the value 
of (Tej / puts a constraint on the range of possible 
values of aH, the semi-hard contribution to the 
hadronic inelastic cross section. The actual mea- 
sured value of (Tej / would give rise to values of (th 
of the order of (Tinei/'^ also at very large cm. en- 
ergies, where one would rather expect an — crinel ■ 
The experimental evidence is also that, in the x 
region accessible experimentally namely at small 
X values, the correlation in fractional momenta is 
not a large effect. 

uh can be worked out rather explicitly when 
only two-body parton correlations are included 
in the many-body parton distributions and when 
each parton can have at most one semi-hard inter- 
action. Two qualitatively different features can 
be present in the two-body parton correlation, 
and both change the relation between gh and 
(Teff with respect to the uncorrelated case: 

1- The distribution in the number of partons is 
not any more a Poissonian, although the depen- 
dence on the kinematical variables of the different 
partons is factorized. 

2- The overall distribution in the number of 
partons, obtained after integrating on the par- 
tonic kinematical variables, is a Poissonian but 
the dependence on the partonic kinematical vari- 
ables is not factorized, in this case the two-body 
parton correlation integrates to zero. 

The general case is obviously a combination of the 
two possibilities. We point out however that both 
cases separately can give rise to a small value of 
(Te// while keeping the value of uh close to 

One can work out explicitly the expression for 
the semi-hard cross section in Eq.(25) considering 
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explicit exainples|ll|. The general result is how- 
ever that the critical value of the impact parame- 
ter Pc, which gives the size to the cross section an , 
is the value which makes small the argument of 
the exponential in the expression of aH{P)- The 
detailed dependence of the argument of the ex- 
ponential at P < Pc is not of great importance 
for the determination of an when, for /3 < /?c, 
the argument of the exponential is already large: 
an is obtained by integrating the probability of 
having at least one semi-hard interaction. When 
the probability to have at least one semi-hard in- 
teraction is close to one, the contribution to the 
integral is very similar for events with the same 
impact parameter and with different but large av- 
erage number of partonic collisions. 

The critical value of the impact parameter 
which gives the size to an is therefore determined 
by the argument of the exponential at the edge 
of the interaction region. The dominant contri- 
bution at the edge is due to the single scattering 
term, since higer order collision terms are impor- 
tant when the density of overlapping matter of 
the two hadrons is large. This is precisely the 
argument of the exponential in the uncorrelated 
case and the consequence is that the resulting 
value of an is not very different with respecy 
to the uncorrelated case. Actually 2ttR^ as dis- 
cussed in the introduction. 

The correlation term is on the contrary able 
to change sizably the effective cross section. One 
may modify the number distribution, without in- 
troducing non-factorized two-body correlations in 
b, by using by using for example the factorized 
expression 

C{u,u') ^ -\D{u)D{u) (30) 
One obtains in this case the relationlllll 



(31) 



If one introduces a correlation term which does 
not modify parton number distribution and which 
therefore integrates to zero, the double scattering 
cross section is incresed, with respect to the un- 
correlated case, by an additive term correspond- 
ing to the convolution of two correlations [O : 



27ri?2 



+ 



+ 



j C(b,b')C(b-/3,b'-/3) 



X (fb(fb'(fp 



(32) 



A qualitative feature is that in both cases one 
obtains a value of cTe// which may be sizably 
smaller with respect to 27ri?2 ~ aH- While, on 
the other hand, nothing prevents the value of an 
from being close to the value of atnei ■ The smaller 
value of CTe/ / , with respect to the expectation of 
the uncorrelated case, is rather generally associ- 
ated with the increased dispersion of the distribu- 
tion in the number of partonic collisions: In the 
case of no correlations the distribution is strictly 
Poissonian when the impact parameter is fixed. 
When correlations are introduced the distribution 
in the number of parton collisions, at fixed /3, is 
not Poissonian any more and the natural conse- 
quence is that the dispersion in the number of 
collisions is increased. 

The indication from the measure of the rate of 
double parton scatterings is therefore that two- 
body parton correlations are likely to be impor- 
tant while, unfortunately, one cannot say much 
about dynamical quantities, like the the correla- 
tion length. Useful observables to be measured, in 
order to get some more insight into the problem, 
would be the semi-hard cross section an and the 
triple parton scattering cross section. The mea- 
sure of an, in association with Ue//, would help 
considerably in clarifying the size of the effect in- 
duced by the presence of the two-body parton cor- 
relations: All present considerations are based on 
the prejudice that an should have a value rather 
close to the value of ainei ■ 

The measure of triple and of higher order par- 
ton scatterings would give important constraints 
on models of the many body parton distributions. 
For example if only lower order correlations where 
important one would be able to fix all the corre- 
lation parameters. 

While a lot of effort has been put in the study 
of the proton structure as a function of the mo- 
mentum fraction x, one should keep in mind that 
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the distribution of partons depends on three de- 
grees of freedom, the momentum fraction x and 
the transverse parton coordinate b. The measure 
of multiple parton collisions is the essential tool 
which allows us to learn on the parton structure 
of the proton in transverse plane. 
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